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Abstract 



Viscoelastic and thermodynamic properties of transient gels comprised of telechelic polymers are theoretically 
studied. We extend classical theories of transient networks so that correlations among polymer chains through the 
network junctions are taken into account. This extension enables us to investigate how rheological quantities, such 
as viscosity and elastic modulus, are affected by the association equilibrium conditions, and how these quantities are 
related to the aggregation number of junctions. We present a theoretical model of transient networks with junctions 
comprised of variable number of hydrophobic groups on the chain ends. Elastically effective chains are defined as the 
chains whose both ends are associated with end groups on other chains. It is shown that the dynamic shear moduli are 
well described in terms of the Maxwell model characterized by a single relaxation time and the high-frequency plateau 
modulus as in the classical theories, but the reduced dynamic shear moduli depend on the polymer concentration and 
temperature through the reduced concentration c given as a combination of the association constant and the volume 
fraction of end groups. The plateau modulus and the zero-shear viscosity rise nonlinearly with increasing c when c is 
small, but they are proportional to c for higher c. The relaxation time also increases as c increases due to the presence 
of pairwise junctions at small c. 
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FIG. 1: (i) Drawing of the "network" postulated in the conventional transient network theories. Each chain interacts with 
fictitious matrix, (ii) Schematic of the network supposed in the present series of papers. Chains interact with each other 
through dissociation/association among end groups. Arrows indicate paths to the (real) network matrix. 



In some polymer gels, junctions can break and recombine in thermal fluctuations or under external forces. They 
are called transient gels or physical gels. Most transient gels exhibit thermoreversible properties, i.e., they reversibly 
change state between gel and sol as thermodynamic conditions vary. Typically, polymers forming such transient 
thermoreversible gels carry a small fraction of interacting groups capable of forming bonds due to associative forces 
such as hydrophobic interaction, ionic association, hydrogen bonding, cross-linking by crystalline segments and so 
on. Among them, hydrophobically-modified water-soluble amphiphilic polymers have attracted widespread interest 
in recent years jlj . Amphiphilic properties stem from the hydrophilicity of the main chain and the hydrophobicity of 
the associative functional groups embedded in the main chain. Attractive force among the functional groups induces 
the formation of transient network in aqueous media above a certain concentration. 

One of the simplest class of associating polymers capable of forming a network is the linear polymers having 
functional groups only at both ends. They are called telechelic polymers. Rheological properties of these polymers 
have been well studied from experimental 0, H, 0, [1, H, [H M, [H, M, M,M,M,m as well as theoretical 
[20l . I21I m, m, m, HE nil point of view with an intention of obtaining a fundamental understanding of associating 
polymer systems. Examples of telechel ic p olymers are poly(ethylene oxide) (PEO) chains end-capped with short 
alkyl groups 0, i, i, i, jl, 0, i, i, M, M,M,M,M,^, perfluoroalkyl end-capped PEO [ul. Hi and telechelic 
poly(A^-isopropylacrylamide) (PNIPAM) carrying octadecyl groups at both ends [18|. They exhibit characteristic 
rheological properties such as temperature- frequency superposition onto a Maxwell fluid 0], breakdown of the Cox- 



In order to investigate molecular origin of these phenomena, Tanaka and Edwards (referred to as TE in the following) 
developed the theory for transient networks [13, EH by extending the kinetic theory for reacting polymers [13] . Under 
the Gaussian chain assumption, TE succeeded to explain, for example, the linear response to the small oscillatory 
deformation described in terms of the Maxwell model with a single relaxation time. Shear thickening can be also 
explained by extending the TE theory so that the tension along the middle chain contains a nonlinear term |26| . We can 
also treat trifunctional associating polymers carrying two different species of functional groups by an straightforward 
extension of the TE theory 28, 29J . Several theories to treat dynamic properties of transient networks have also been 
presented up to now. For example, Wang [12] took isolated chains into consideration, and Vaccaro and Marrucci [13] 
incorporated the effect of incomplete relaxation of detached chains. 

In all transient network theories proposed so far [H, Hlj, [13, [H, [13j [H) [H, [1^ , it is implicitly assumed that a 
fictitious network exists a priori (see Fig|T] (i)). This network matrix is not a substantial one in a sense that it 
itself does not contribute to the elasticity of the system, but it plays a role as a substrate of the chains on which 
association/dissociation of the end groups takes place. Chains whose both ends are connected with this matrix 
are supposed to be elastically effective. Correlations among chains are not taken into account in this treatment 
because each chain interacts only with such matrix, and consequently the concentration dependence of the rheological 
quantities cannot be properly predicted; the elastic modulus, viscosity etc are simply proportional to the polymer 
concentration. Furthermore, it is difficult by definition to incorporate the information about the network junction 
such as the aggregation number and to study the effect of surfactants added to the system. 

In this series of papers, we remove this unfavorable assumption and develop a theory of thermoreversible transient 
networks formed by multiple junctions comprised of any number of functional groups as depicted in Fig[T](ii). This 
modification enable us to deal with the sol/gel transition phenomenon and to predict the proper concentration 
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dependence of the dynamic shear moduh G' and G" in the postgel regime of the solution. In the following, we refer to 
the number of functional groups per a junction (so called aggregation number) as the junction multiplicity according 
toref. [13]. 

In this paper, we formulate a theory of transient network formed by thermoreversiblc multiple junction. Elastically 
effective chains, or active chains, are defined only in local, i.e., chains whose both ends are associated with other 
functional groups are supposed to be elastically effective. We derive a formula to calculate G" and G" as a function 
of the frequency w of applied small oscillation, and analyze how G" and G" characterized by the high-frequency 
plateau modulus Goo, zero-shear viscosity rjQ and relaxation time r depend on the polymer concentration. It is 
shown that Goo and rjo increase nonlinearly with concentration at lower concentration and t also increases with 
concentration corresponding to experimental observation [1, [2^ . In equilibrium, we find the following mass action 
law: tpk = K]^{il)i)^ under the assumption that the connection rate of a functional group to a junction with the 
multiplicity k (called /c-junction) is proportional to the number of functional groups in /c-junctions. {tpk is the volume 
fraction of the functional groups belonging to fc-junctions, and Kk is the reaction constant for the formation of a 
A:-junction from k isolated functional groups.) This relation is equivalent to the multiple-equilibrium condition that 
Tanaka and Stockmayer derived in the theory of associating polymer solutions [30| . In order to investigate the critical 
behavior of G' and G" near the gelation point, we need to define active chains taking account of global structure of 
the network. Studies on the basis of this definition for active chains will be presented in the following paper. Looped 
chains, and then flower micelles comprised of loops, are assumed to be absent for simplicity throughout this series. 
Effects of the competition between the intra- and intermolecular association at a junction as well as the effect of 
additives such as surfactants [Sll, [H, [sl] will be investigated in the forthcoming papers. 

This paper is organized as follows. In section [TTl we will derive time development equations for the distribution 
function of chains whose one end is incorporated into a fc-junction while the other end is belonging to a fc'-junction. 
Kinetic equations for these chains will be also derived in this section. In section IIIIl we will introduce the dissocia- 
tion/association rates of functional groups. Equilibrium distribution of the junction multiplicity and its mean value 
as a function of polymer concentration will be discussed in section IIVI In section |Vl we will study linear rheology of 
the system. Finally, a relation between the present theory and the TE theory will be discussed in section IVll 

II. TIME DEVELOPMENT OF TRANSIENT NETWORKS FORMED BY JUNCTIONS WITH 

VARIABLE MULTIPLICITY 

A. Assumptions 

We consider athermal solutions of linear polymers (or primary chains) carrying two functional groups at their 
both ends. Here, functional group is a groups or a short segment of the primary chain that can form aggregates (or 
junctions) in a solution through the noncovalent bonding. Primary chains can associate with each other in a solution 
through the aggregation of functional groups, while they can be detached from others due to thermal agitation or 
macroscopic deformations applied to the system. In equilibrium, thermodynamic conditions (temperature, polymer 
concentration etc) determine the association/dissociation rate of functional groups and hence the amount of junctions 
in the solution. In this sense, these junctions are thermoreversible. Above a certain polymer concentration, primary 
chains eventually construct a macroscopic network physically cross-linked by these junctions. Thus the network we 
are going to study is also thermoreversible as well as transient. The association-dissociation reaction of functional 
groups is assumed to occur step by step. 

We allow junctions to be formed by any number of functional groups. According to Tanaka and Stockmayer [30| . 
the number of functional groups forming a junction is referred to as the junction multiplicity. We also call the junction 
of the multiplicity k (= 1, 2, 3, • • • ) the fc-junction, i.e., 1-junction is an unreacted group (or a functional group which 
remains unassociated) , 2-junction is a pairwisc junction, etc. For a while, let us identify a head and a tail of each chain 
for convenience by marking one end of each chain temporarily. Of course, this does not affect mechanical properties of 
our system. Then we term the primary chain whose head is incorporated into a fc-junction whereas tail is a member 
of the distinct fc'-junction as the (fc, fc')-chain. For instance, (fc, l)-chain (fc>2) is a (primary) dangling chain whose 
tail is not connected with any other groups, and (1, l)-chain is an isolated chain. 

We assume that chains are elastically effective when their both ends are connected with any other functional groups. 
These chains are called active chains. Active chains are assumed to deform affinely to the macroscopic deformation 
applied to the system. Note that active chains are defined only in local in a sense that chains are elastically effective 
irrespective of whether functional groups they are connected with belong to the infinite network or not. As shown 
below, we can study the sol/gel transition phenomenon only in unsatisfactory in this theoretical framework. The 
definition of active chains will be extended in the following paper so that global information of the network can be 
included. The Rouse relaxation time of the primary chain is assumed to be much smaller than the characteristic 
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time of a macroscopic deformation applied to the system and the hfetime of active chains, so that chains in elastically 
non-effective (i.e., dangUng and isolated) states are virtually in an equilibrium state even under flow. The primary 
chains are assumed to be Gaussian with uniform molecular weight M (number of repeat units is N) that is smaller 
than the entanglement molecular weight. 



B. Time-Development Equation for Active Chains 



Let Fk^k' (r, t)dT be the number of {k, fc')-chains at time t per unit volume having the head-to-tail vector r r + dr. 
Then the total number Vk,k'{t) of {k, fc')-chains (per unit volume) is given by 

i^k,k'{t)^ Jdv Fk,k'{r,t). (II.l) 

Dangling and isolated chains are the Gaussian chain irrespective of time, that is, 

Fi^k'{r,t) = iyi,k'it)foir) (for fc > 1), (II.2a) 

Fk,iir,t)^iykAt)foir) (for A:' > 1), (n.2b) 



where 

3/2 

2TTNa^ ) ^''"^ V 2A^a2 J 

is the probability distribution function (PDF) that the chain takes the end-to-end vector r (a is the length of a repeat 
unit of the chain). The number of chains whose head is incorporated into a fc-junction is given by 



/o(r) ^ ( 7^ ) exp ( (II.3) 



x1"W = E^mW, (11.4) 



v-; — 
i>i 

while the number of chains whose tail is belonging to a fc-junction is 

1>1 

Then the number of chains whose one end, irrespective of whether it is a head or tail, is incorporated into a fc-junction 
(called fc-chain hereafter) can be expressed as 

Xk{t) = l{x^'\t)+x'-k\t)). (II.6) 

A factor 1/2 is necessary to avoid double counting of the fc-chain. Our aim in this section is to find time-development 

equations that Fk^k'i^it), I'k.k'i'^) ^^^d X^'*(t) (dot represents h or t) obey when a macroscopic deformation represented 
by the rate of deformation tensor k{t) is applied to the system. 

Firstly, we derive the equation for active (fc, fc')-chains (fc, fc' > 2). According to the affine deformation assumption, 
Ffe^fc/(r, i) satisfies the following evolution equation: 

dFk,k^ij, + V ■ {v{t)Fk,k' (r, t)) ^ Wl!:l (r, t) + wi%, (r, t) (forfc,fc'>2), (II.7) 

where v{t) = k{t)r is the rate of deformation of the head-to-tail vector r of the active chain. The first (or the second) 
term in the right-hand side of (|II.7|) is a reaction term that describes the net increment of Fk^k' (r, t) per unit time 
due to association-dissociation reactions among the head (or tail) of the (fc, fc')-chain and the functional groups on 
the other chains. The kinetic equation for active (fc, fc')-chains is then given by integrating both sides of pi.7p with 
respect to r. We find 

where 

4 ) {t) ^ [dr Wl% (r, t) (for fc, fc' > 2). (II.9) 
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FIG. 2: Reaction between the functional group on the {k, fe')-chain and the functional group on the other chains. Circles 
indicate junctions and lines emanating from circles represent the primary chains. The letter inside the circle denotes the 
junction multiplicity. Smaller circles with the number 1 inside represent unreacted functional groups. Chains depicted by bold 
lines participate in the reaction. The number of (fc, fc')-chains change if one of these (stepwise) reaction occur. Reaction rates 
are shown near arrows for each reaction. 



The reaction terms M^^^, (r,t) and w^'^, (t) are derived according to the following procedure. The number of active 
{k, fc')-chains decreases if** 

(i) the head of the {k, fc')-chain is dissociated from a fc-junction, or 

(ii) a functional group on the other chain is dissociated from the head of the (k, fc')-chain, or 

(iii) an unreacted functional group connects with the head of the (fc, A:')-chain. 
On the other hand, the number of active (fc, A:')-chains increases if 

(iv) the unreacted head of the (l,fc')-chain is connected with fc — 1-junction, or 

(v) the unreacted group of the other chain is connected with a head of the {k—1, fc')-chain, or 

(vi) a functional group on the other chain is disconnected from the head of the (fc+1, /c')-chain. 

These reactions (i) ~ (vi) are schematically depicted in FigO where the corresponding reactions regarding the tail of 
the (fc, fc')-chain, (i') ~ (vi'), are also shown. Now we introduce the probability (or dissociation rate) Pk{r) that a 
functional group incorporated into the fc-junction (fc > 2) detaches itself from the junction per unit time. In general, 
it is an increasing function with respect to the end-to-end length r ~ |r| of the fc-chain [20|. Then the increment of 
(fc, fc')-chains with the head-to-tail vector r due to the reactions (i) and (ii) per unit time is written as 

•' s=l s' = l \ s" = l s" = l / 

where the first term Pki^) in the parenthesis stems from the reaction (i) while the second and the third terms are 
originated from the reaction (ii). Tk^k' (r, {^^'^^, {r^*^}; t) is the number of (fc, fc')-chains with the head-to-tail vector r 
whose head is incorporated into the fc-junction formed by kh heads of the other chains each having the head-to- vector 
r^'''', ■ • ■ , r^^^ and kt tails of the other chains each having the head-to- vector r^*\ ■ • ■ ,r^i^^ (see Figl3]). Note that a 
relation {kh + 1) + kt — k holds. We can express J-'k,k' as follows: 

J-fc,fc,(r ,{r('')},{rl*)};t) = Cfe,i^,,,-(r,t)n/f)(rW,t) f[ fi'\r^^ ,t), (ILll) 

s=l s' = l 



We are considering reactions regarding the head of the (fc, fe')-chain. The multipHcity of the junction to which its tail is connected is 
assumed to be fixed at fc'. 
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:^,A'(r,{rW},{r«}) 




FIG. 3: Explanation of eq. HII.10|I . Circles stand for junctions. Arrows connecting two junctions represent the head-to-tail 
vector of the chain (chains are not shown for simplicity). The number of the state A is given by The first term in 

the parenthesis of (|n.lOp describes the transition from the state A to the state B, while the second (or the third) term in the 
parenthesis of (|II.10|I represents the transition from A to C (or to D). 



where 



j:i>iFkAr,t) 



(11.12) 



is the PDF that the chain whose head is incorporated into a /c-junction has the head-to-tail vector r, and similarly 



Jk y'^^'^) - - rpr- 



(11.13) 



is the PDF that the chain whose tail is belonging to a fc-junction takes the head-to-tail vector r. Note that ■'(r, t) 
is normalized to one, i.e., 



Jdrfl:\r,t) 



1. 



(11.14) 



A prefactor Ck^ is the probability that the fc-junction is formed by kh heads and kt tails (in addition to a head of the 
(fc, fc')-chain), and hence it is given by 



1 (fc-1)! 



2^-1 kh\kt\ ■ 

Upon substitution of prTTjl together with (|ILT2|1 . (|ILT3)) and (|ILT5)) . eq. (|irTO)) reduces to 

Pk{r)Fk,k'{r,t) + {k- l){f]k{r)){t)Fk,k' (r,t), 

where 



mr)){t) = J dvh{r)h{v,t) 
is the expectation value of /3fc (r) averaged with respect to the PDF for fc-chains given by 

/,(r,t)^i(/f)(r,i)+/f (r,t)). 



(11.15) 



(11.16) 



(11.17) 



(11.18) 
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The number of {k, fc')-chains decreases through the reaction (iii) only when the fc-junction, to which an unreacted 
group are going to connect, includes the head of a (fc,fc')-chain (with head-to- vector r). Here let us introduce the 
probability (or connection rate) Pk{t) that an unreacted functional group catches a fc-junction per unit time.*^ Then 

the decrement of Fk,k'{Y,t) due to the reaction (iii) is written as the product of the number Pk{Xi\^) + Xi\^)) of 
unreacted ends (both head and tail) that catches the fc-junction per unit time and the number of (fc, fc')-chain with 
head-to-vector r per a fc-junction, that is, 

Pkit) (x^t\t)+xf{t)) ^ kpk{t)^^Fk^k'{v,t), (11.19) 

where 

is the number of fc-junctions and a relation (|II.6p has been used for the equality. The increment of the (fc, fc')-chain 
with the head-to-tail vector r through the reaction (iv) is given by 

Pk-iit)Fi,k'ir,t). (11.21) 

The increment of Fk^k' (r, t) due to the reaction (v) is obtained by replacing fc in pi.l9p with fc — 1 as follows: 

{k~l)pk-i{t)^^Fk-i.k'{r,t), (11.22) 
Xk-i{t) 

whereas the increment due to the reaction (vi) is given by replacing fc in the second term of pi.l6p with fc + 1, i.e., 

k{pk+i{r)){t)Fk+i,k'ir,t). (11.23) 
Summarizing, the reaction term regarding the head of the (fc, fc')-chains found to be 

W'S'Ci",^ =-Pk{r)Fk^k'{r,t) +pk-i{t)Fi^k'{r,t) 

~Bk{t)Fk,k' (r, t) + Bk+i {t)Fk+i,k' (r, t) 

-Pk{t)Fk,k'{r,t)+Pk^i{t)Fk^i,k'{r,t) (for fc,fc' > 2) (11.24) 



where we have put 



Bk{t) = (fc-l)(/3fe(r))(i), (IL25a) 

Pkit) = kpkit)^. (II.25b) 
Xk{t) 



Integrating (|II.24p with respect to r, we also obtain'^ 



w 



CO 

k,k' 



(t) =- / dr /3k{r)Fk,k'{r,t) + Pk-i(t)vi,k'{t) 



~Bk{t)vk,k'{t) + Bk+i{t)vk+i,k'{t) 

~Pk{t)vk,k'{t) + Pk-i{t)vk^i,k'{t) (for fc > 2). (11.26) 
According to the similar procedure, reaction terms regarding the tail of the (fc, fc')-chain are derived as follows: 

Wl!l.{v,t) ^-(3k'{r)FkM'{v,t) +pk'-i{t)Fk,i{v,t) 

-Bk' {t)Fk,k' (r, t) + Bk.+i{t)Fk,k'+i (r, t) 

-Pk'{t)Fk,k'{r,t) + Pk--i{t)Fk,k'-i{r,t) (forfc,fc'>2) (11.27) 



^ In general, Pk{t) depends on the number /ifc(t) of the fc-junction in the system; the more the fc-junction, the higher the probabiUty for 
an unreacted group to catch the junction (see section [llll for detail). Thus Pk{t) depends on time through fik{t)- 
Eq. llTL26t holds for fc' > 1, while l llL28l l holds for fc > 1. 
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and 



-Bk'{t)vk,k'{t) + Bk'+i{t)vk,k'+i{t) 

-Pk'{t)vk^k'{t)+Pk'-i{t)uk,k'-i{t) (forfc'>2). (11.28) 

Next, we derive the kinetic equation for dangling chains, i.e., eq. pi.8[) with fc = 1, fc' > 2 or fc > 2, fc' = 1. The 
reaction term wj'*^, (t) can be derived as foUows. The (1, fc')-chain is generated by detaching a head of the (Z, fc')-chain 
{I > 3) from the Z-junction. The increment of i>i_k' due to this process per unit time is 

J2 [ dv Pi{r)Fi^k'{v,t). (11.29) 

The number of the (1, fc')-chain also increases when the 2-junction to which a head of the (2, /c')-chain is belonging 
divides into two unreacted groups. The increment according to this procedure is given by (see Appendix E]) 

dv f32{r)F2,k'{r,t) + (^J dr I32{r)h{^,t)^ V2,k'{t). (11.30) 

On the contrary, (1, fc')-chains are annihilated if the head of the (1, /c')-chain is connected with Z-junction {I > 2). The 

decrement of vi^k'it) caused by this process is given by (^i>2Pii^)j^^,k'it)- Besides, (1, A;')-chains are annihilated 

when its head captures the unreacted group. The decrement due to this reaction is estimated to be 2pi{t)i'i^k'it) (see 
FigHllin Appendix lAl as a reference) As a result, we obtain the reaction term with respect to the head of (1, A:')-chain 



as 



l>2 " 1>1 

+ B2{t)v2^k'(t) - pi{t)vi^k'{t). (11.31) 

The reaction term u}-^\, regarding the tail of (1, fc')-chain (fc' > 2) is given by (|II.28P with k ~ 1. According to the 
similar way, the reaction term for the tail of (fc, l)-chains is found to be 



(t) 

w 



t>2 " l>l 

+B2{t)vk,2{t) - pi{t)vk,i{t) (11.32) 

and (k > 2) is given by (jll.28|l with k' — 1. Finally, the kinetic equation for isolated chain can be obtained by 
setting fc = fc' = 1 in jlLSl) . where w^'i (i) (or wj*l(i)) is given by (|lL3l|) with k' = I (or k = 1). 

By making use of (|!L26|1 . (HlMJ, (|IL3T|) and (|lL32l) . we can also obtain the equation for xVi^) follows: 

(^1 dr (3kir)fi-\r,t)^ xi\t) + pk-^it)x'i\t) 
-Bk{t)xi\t) + Bk+i{t)xl:Ut) 

-Pkit)xi\t) + Pk-i{t)xil,{t) (for fc > 2), (II.33a) 

l>2 ' 1>1 

+B2it)x^2\t)-Pi(t)x[\t). (n.33b) 



Eq. lHOTI l holds for any k'{> 1) and l llL32l l is satisfied for any A;(> 1). 
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One can confirm from pi.33|) that the total number of chains conserves, i.e., 

fc>lfc'>l fe>l k'>l 

Note that pi.34p holds for arbitrary /3fc(r) and pk- Hereafter, we denote the number of total chains (per unit volume) 
as n, i.e., 



fe>i fe'>i 

Up to now, the head and the tail of the chain have been distinguished for convenience. Since we are treating 
symmetric chains actually, however, subscript of h'k,k' is exchangeable: i'k,k'{t) ~ i^k,k'it), and hence kinetic equations 
for (fc, A:')-chain (|IL8)) with the reaction term (|IL28|) . 0^ and ' (|I02)) can be summarized into 

^I^bK^ = Wk.k'{t)+Wk-Mt) (11.36) 
at 

where, for fc' > 1, 



Wk,k'it) = - j dr(3kir)Fk,k'{r,t)+pk-i{t)i^iM'{t) - {Bk{t) + Pk{t)W,k'{t) 

+ Bk+i{t)vk+i,k'(t)+Pk-i{t)vk-iM'{t) (for fc>2), (II.37a) 



.k'{t)^Y. I dv Pi{r)Fuk'{v,t) + B2{t)v2,k'{t) ~ Pi(0 + E^^'(^) ^^,k'{t). (II.37b) 

i>2 \ '>1 / 



Further, we find that Xk{t) — Xk^H^) — Xk H^) satisfies the following equation: 

dXkit) 



dt 

where 



ukit), (11.38) 



Ukit) =-k{f3k{r)){t)xk{t) + K(ik+i{r)){t)xk+i{t) 

+kpk-i{t)xi{t) - kpk{t)xi{t) (for k > 2), (II.39a) 

Mt) =^(A(r))(t)xz(i) - (j2P^it))xi{t) 



l>2 1>1 

+ (/32(r))(i)x2(i) -Piit)xiit). (n.39b) 

Evolution equation pi.7[) with the reaction terms given by pi.24[) and pi.27[) as well as the kinetic equations for 
(fc, fc')-chain pi.36|) together with pi.37p and for fc-chain pi.38p together with pi.39|) are fundamental equations of 
our transient networks and will be solved in the following sections. 



III. REACTION RATES OF FUNCTIONAL GROUPS 



In the rest of this article, we treat the dissociation rate of the functional group as constant independent of the 
end-to-end length of a chain. This treatment is valid since we are going to apply a small oscillatory deformation to 
the system and the change of the dissociation rate through r is quite small in this case. Then, according to TE, 
the dissociation rate from the fc-junction takes a form: f3k = exp (— VFfc/fcsT) [2l| where luq is a reciprocal of a 
microscopic time and Wk is a potential barrier for the dissociation. We further assume that the potential barrier does 
not depend on the multiplicity of the junction and set Wk = W for all k. Then the dissociation rate also does not 
depend on the junction multiplicity and is written as 



Pk^ujf^expi-W/kBT) (=/?). 



(III.1) 
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FIG. 4; The potential barrier around a fc-junction for the dissociation-association reaction of the functional group. 



The connection rate of an unreacted group to a fc-junction should increase with increasing the number of functional 
groups forming the fc-junction in the close vicinity of the unreacted group. We assume that the connection rate takes 
the form 



Pk{t) = ujo exp {-{W - e)/kBT) k^ik{t)vohk, (III.2) 

where e is a binding energy between the functional group and the junction'^ (see Fig|4]), kfik{t)vo is the number of 
functional groups forming fc-junctions in the effective volume vo of the (unreacted) functional group, and hk is a 
proportional factor given in the next section. Note that the connection rate depends on time through the number of 
fc-junction in general. Eq. pil.2p can be rewritten as 

Pkit) ^ f3X{T)4,qk{t)hk (III.3) 

where 

is the probability that an arbitrary chosen functional group belongs to a fc-junction, ip = 2nvo is the volume fraction 
of functional groups, and 

A(T) ^ exp (^^) (in.5) 

is the association constant of the reactions introduced in ref. [30l |. Thus the connection rate to the fc-junction is 
proportional to the volume fraction ipq^. of fc-junctions. 



IV. EQUILIBRIUM PROPERTIES 



In equilibrium state, the number Xk of fc-chain, or equivalently qk, is formally obtained by setting the left-hand side 
of pi.33[) to 0. (Here and hereafter, all quantities in equilibrium are denoted without the argument t. For example, 
Vkjk' is the number of (fc, fc')-chains in equilibrium.) Thus we find 

qk^jfj^qi (forfc>2) (IV.l) 



where qi is obtained from the normalization condition X]fc>i 9fe = 1 



,f 



^ + 2lk>2Pk-l/{Pk{r)) 



It is also assumed that the binding energy does not depend on the junction multiplicity. 
^ This equation is equivalent to the number conservation equation X)fe>i Xk = " for chains. 
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On substitution of (IIII.ip and pil.3|) . eq. (jlV.ip becomes 

qk = Xijjhk-iqk-iqi (for k > 2). (IV.3) 

By an iterating procedure, we find the following mass action law: 

gfc -7fe(AV)'"'gf (forfc>2), (IV.4) 

qi - l/liz), (IV.5) 

where 

7(z)^^7fc^'-' (IV.7) 

k>l 

and z = \ipqi . If A and ip are given, gi can be derived by solving (jIV.SP . Then we can obtain qk [k > 2) from (jrO]) . 
The association condition (IIV.4P together with (IIV.SP have been derived in the theory of thermoreversible gelation 
with junctions of variable multiplicity developed by Tanaka and Stockmayer (referred to as TS) [30|.^ In the TS 
theory, jk is interpreted as a factor giving the surface correction for the binding energy, although it is set to unity 
for all k for simphcity. We will adjust hk (and hence 7fc) to derive specific models for junctions (see below). TS has 
shown that most quantities describing transient gels in equilibrium depend on the polymer volume fraction ■0 through 
the combination of A(r) and tp{= 2(f)/ N). This holds not only in equilibrium state but also under small deformation 
as shown in the next section. Therefore, we use c = X(T)'ip as the reduced polymer concentration in the following. In 
equilibrium, ()II.37P reduces to a simpler form: 

Wfe.fc' = -/3fc(l + zhk)yk,k' + Pkuk+i,k' + P{k - l)zhk-iiyk~i,k' + f3chk^iqk-ivi,k' 

(for fc > 2), (IV.8a) 

= /3 J/2,fc' + ^ t'/^fc' - I3c{hiqi+^^hiqi\vi,k'. (IV.8b) 

y i>2 J ^ i>i ' 

By solving a equation {dvk^k' /dt=) Wk,k' + uik\k = 0, we can obtain the number of (/c, fc')-chains in equilibrium as 

i^k.k' = nqkqk'- (IV.9) 

Note that Vk.k' /n = qkqk' is the probability that one end of an arbitrary chosen chain is incorporated into a fc-junction 
while its other end is belonging to a /c'-junction. Thus (jIV.9[) is compatible with the definition of the (fc, fc')-chain. 

Now we consider two special cases as for the multiplicity that the junction can take; 1) a saturating junction model 
and 2) a fixed multiplicity model. These two models have been considered by Tanaka and Stockmayer [sO] in studies 
of the phase behavior of associating polymer solutions in equilibrium.^ In the saturating junction model, junction 
multiplicity has an upper limit s™, that is, the multiplicity is allowed to take a limited range fc = 1, 2, ...Sm- This model 
might be applied to the junction formed by ionic dipolar interaction or hydrophobic aggregation because the space 
around such a junction is packed with polymer chains and steric hindrance among them restricts the number of chains 
connected with the junction. [30| . The mean multiplicity generally depends on the reduced polymer concentration 
in this model. On the other hand, in the fixed multiplicity model, the multiplicity can take only one fixed number s, 
i.e., we have only fc = 1 (unreacted) and k = s (reacted) irrespective of the reduced polymer concentration. 

A. Saturating Junction Model 

We can impose the upper limit on the junction multiplicity by assuming that hk is given as a stepwise function: 

1 (l|fc<.„-l) ^ 
(fc > s„) 



^ In this series of articles, we denote the probability that an arbitrary chosen functional group belongs to a fe-junction as gj,, although the 

same quantity is represented as pf, in ref. '30| , because the symbol pj, is used as the connection rate of functional groups. 
^ They have also considered a minimum junction model [3Cll |. where junctions are allowed to take the multiplicity k = I (unreacted) and 
= so, so + 1, ■ • • . 
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5^-4 5,„-12 5m = 20 




FIG. 5: The probability distribution qt that a randomly selected functional group to be in a fc-junction for the maximum 
multiplicity Sm = 4 (left column), 12 (middle column), 20 (right column), and for the reduced polymer concentration c — 1 
(top raw), c — s,n (middle row), c — 1000 (bottom row). The weight-average multiplicity of the junction is shown in each 
figure. 



For such hk, PV.6|) reduces to 



and hence qk takes the form 



Jl (l<fc<s™) 

^'^^jo (fc>.™ + i) ' (^^-^^^ 



Thus, junctions with the multiphcity larger than s,,,, do not exist anymore. We can obtain qi by solving pV.5|) : 

1 _ 1 - (cgi)^'" 



qi 1 - cqi 



(IV. 13) 



Note that the right-hand side of (|IV.13[) (denoted as g{qi) for simplicity) is Sm when qi — 1/c. Therefore, when the 
condition = c is satisfied, a solution of (|IV.13[) is qi = 1/c. As a result, qk does not depend on k, i.e., q^- = 1/c 
for all fc(< Sm) (see middle row figures of Fig[5]). When Sm is larger than c, g{qi = 1/c) = Sm is also larger than c. 
Therefore, a solution of (jIV.13|) satisfies a condition gi < 1/c because g{qi) is an increasing function with respect to 
qi. Consequently, qk is a decreasing function with respect to k (see top raw figures of FigU]). On the contrary, when 
Sm is smaller than c, a solution of (jIV.13p fulfills a condition qi > 1/c, and hence qk is an increasing function of k 
(see bottom row figures of FiglS]). 

FigUl (i) shows the extent of reaction a = 1 — gi as a function of for several reduced concentration c. We see 
that a approaches to a fixed value for each c as Sm increases. This value can be estimated as follows. In an extreme 
case where Sm is much larger than c, the right-hand side of pV.lSp is close to 1/(1 — cqi) thanks to the condition 
cqi < 1. Thus we find 

91 = -^ (forSm»c), (IV. 14) 



and therefore 



In this extreme, qk is expressed as 



1 + c 



(for Sm > c). (IV. 15) 



(IV. 16) 



13 



(i) (ii) 




10 20 30 40 5 10 15 20 

Sff, c 



FIG. 6: The extent of reaction a of the saturating junction model as a function of the maximum multiplicity Sm (i), and of 
the reduced polymer concentration c (ii). The reduced concentration is varying from curve to curve in (i), while the maximum 
multiplicity is changing in (ii). The inset of (ii) shows the linear- log plot of a as a function of c. Dotted curves (behind the 
curves for Sm ~ 20) in (ii) represent a = c/{l + c). 



(i) (ii) 




FIG. 7: The weight- average multiplicity fiw of the junction for the saturating junction model as a function of the maximum 
multiplicity Sm (i), and of the reduced polymer concentration c (ii). The reduced concentration is varying from curve to curve 
in (i), while the maximum multiplicity is changing in (ii). The inset of (ii) shows the linear- log plot of fiw as a function of c. 



where 

K = l/log[(l + c)/c] (IV. 17) 

indicates the width of the distribution. As an example, (|IV.16|) is drawn in FiglS] for c = 1 and Sm — 20. As seen 
from Fig|S] (ii), the extent of reaction behaves as c/(l + c) for c much smaller than Sm, while it approaches to 1 if c 
exceeds Sm- 

The weight-average multiplicity defined by 

M™ = £fcgfc (IV. 18) 

fc=i 

is shown in Fig|5]for each s„i and c and is depicted in Fig|7]as a function of Sm (i) and of c (ii). Note that /i^ includes 
unreacted groups as fc = 1 junctions. When Sm is much larger than c, it is close to (see FiglT] (i) and (ii)) 

/i^ = 1 + c (for s,n > c). (IV. 19) 

When the reduced concentration is low and satisfies the condition c ^ we see from top three figures of Fig[5]that 
fiw is much smaller than Sm and close to 1 irrespective the value of s,„. This is because there are many unreacted 
groups (represented by qi) in the system. On the other hand, when the reduced concentration is high and c ^ Sm, 
is close to Sm (see bottom three figures of Fig[5]) since many functional groups are incorporated into junctions 
with k = Sm- These tendencies can also be confirmed from FiglT] (i). When Sm is much smaller than c, then /Xu, is 
approximately Sm, while /x^, approaches to 1 -I- c when Sm exceeds c for each c. As seen from FiglT] (ii), /iu, behaves 
as 1 -I- c irrespective of the value of Sm when c <C Sm, and it approaches to Sm as c increases beyond Sm- 
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FIG. 8: The probability distribution qk that a randomly selected functional group to be in a fc-junction for the reduced 
concentration c = 1 (left), 2 (middle), and 10 (right). The multiplicity of the junction is fixed at s = 10. 



B. Fixed Multiplicity Model 

Let hk be 



5 (1 < fc < s - 1) 

hk = { 5^^'-^^ {k^s-1) , (IV.20) 



{k>s-l) 



where 5 is a positive value. Then qk becomes 



iSc)''-^q'l (1 < fc < s) 
qk = { c'-'qf (fc = s) . (IV.21) 

(fc > s) 

All junctions take approximately the same multiplicity s if we employ 6 much smaller than 1' because qk is nearly 
equal to except for the case of fc = s (and fc = 1), i.e., 

Qk^f^f'^' [l^l] (for<5«l). (IV.22) 

It should be emphasized here that we cannot fbc the junction multiplicity rigorously at s by putting 6 = because 
junctions whose multiplicity is less than s must exist for the creation of s-junctions under the assumption of stepwise 
reactions. We set 6 = 0.01 in this series of papers. This i5 value is enough to describe characteristic properties of the 
system where junctions can take only one multiplicity. In the following, we often use the equal sign instead of the 
nearly equal sign (~) for equations (such as pV.22[) ) that approximately hold for small S. In the fixed multiplicity 
model, the extent of reaction is given hy a = qs because of the normalization condition qi + qs = I. 
The probability qi to find an unreacted group can be obtained by solving (IIV.SP : 

- = l + {cq,r~\ (IV.23) 

The right-hand side of PV.23P (denoted as g{qi)) is equal to 2 for qi = 1/c. Therefore, if c = 2, a solution of 
PV.23P is given by qi = l/2(= q^). If c > 2, a solution of pV.23p satisfies a condition 1/c < qi < 1/2 because 
g{qi = 1/2) = 1 + {c/2)^^~^'> is greater than 2. Then qs{> 1/2) is larger than qi (see FigE]) indicating that there are 
more reacted groups than unreacted ones. Note that when s is much larger than 1, g{qi) is equal to 1 for qi < 1/c but 
it almost diverges at gi > 1/c, and therefore, a solution of PV.23|) approaches to 1/c. On the other hand, if c < 2, the 
solution of pV.23p satisfies a condition 1/2 < gi < min(l, 1/c), and it approaches to 1/c for c satisfying 1 < c < 2, or 
to 1 for c < 1 with increasing s. Summarizing, 

'?i = {l^' (c<l) (1^-24) 



or, equivalently, 



I > Jj (for . » 1). (IV.25) 



' On the other hand, the saturating junction model can be realized by setting 5 = 1 in I IIV.201 1. 
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FIG. 9: The extent of reaction of the fixed multiplicity model as a function of the reduced polymer concentration c (i), and of 
the junction multiplicity s (ii). The inset of (ii) shows the linear-log plot of a as a function of c. 



Thus junctions suddenly appear at c = 1 in this extreme. 

Fig[9] (i) shows the extent of reaction as a function of s. We can confirm the tendency stated above. If c < 1, the 
extent of reaction approaches to with increasing s. If c > 1, on the other hand, it approaches to 1 — 1/c for a given 
c. But when c = 2, it does not depend on s. FiglQ] (ii) shows the extent of reaction as a function of c. It increases 
abruptly around c — 1 when s is large. Such a sharp increase in a stems from the fact that the junctions can take 
(approximately) only one multiplicity; even if several functional groups spend a certain duration of time in the close 
vicinity of each other, they cannot aggregate unless s groups participate in this event. Actually, this tendency cannot 
be observed in the saturating junction model (see Fig|S]) since the junction can be formed by any number of functional 
groups less than a certain value. 



V. DYNAMIC-MECHANICAL AND VISCOELASTIC PROPERTIES 



Now, we apply a small oscillatory shear deformation to the present system whose rate of deformation tensor is 
represented by 

(0 euj cosmt \ 
, (V.l) 
y 

where e is a dimensionless infinitesimal amplitude and uj is the frequency of the oscillation. By substituting (jV.ip into 
the left-hand side of (|II.7p . and by putting (jllLip and (|IIL3p into the right-hand side (given by pi.24p and pi. 271) ). 
we have the evolution equation for Fk^k' (r, t) given as follows: 

dFk,k'{r,t) dFk,k'{r,t) 

dt + dx ^^"^"^"^ 

= -Qk,k'{t)Fk,k'{v,t) + Bk+iFk+iM'{'c,t) + Bk'+iFkM'+ii'c.t) 

+Pk-i{t)Fk-i^k'{v, t)+Pk'-i{t)Fk,k'-i{v,t) 

+pc{hk-iqk-iit)vk',iit) + hk'-iPk'-i{t)vk,i{t))foir) (for k,k' > 2), (V.2) 



where^ 



as defined in (jll.25p . and 



Bk = Pik-l), (V.3) 
Pk{t) = I3z{t)khk (V.4) 



Qk,k'{t)=Pk + Pk{t)+pk' + Pk'{t). (V.5) 



j z{t) = \iPqi(t). 
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Here, we expand F^^k' (r, t) with respect to the powers of e up to the first order: 

Fk,k' (r, t) = Fjfl (r) + eFj^'l (r, t). (V.6) 

On the other hand, the first-order terms of i^k,k'{t) and qk{t) are neghgible [IOjUIIj i-e-, they have equihbrium values 
and satisfy the association condition (jIV.4p even under the small shear deformation. The zeroth-order term of Fk,k' (r) 
represents its equilibrium value in the absence of the oscillation, and hence it is written as 

Fi%ir)^i.k,k'Mr). (V.7) 
By comparing the order, we see that the first-order term Fj,^l,{r, t) satisfies the following equation: 

m i^k.k'j^foirjujcosujt 

= -Qk,k'F^%{r,t)+Bk+,Ft';>,,,{r,t)+Bk-+iFt%^,{r,t) 

+Pk-iFi%^^,{v,t)+Pk'-iFl^l,_^{Y,t) (forfc,A:'>2). (V.8) 



k 



The solution of (|V.8[) after the long-time limit takes the form' 

= (ff^,Msina;t-f 5l'^,,Hcosc^t)-^/o(r). (V.9) 

Note that subscripts of g' (and g") are exchangeable, i.e., g'j, j,, = g'j,, f. (or g'l ^ = g'y j.), and we can put g'^. ^~g'l^=Q 
for fc > 1. Substituting (|V.9p into (jV.8[) . we obtain the simultaneous equation for g'^. f., and g'^^., as follows: 

9'k,k' - {-QkM'glk' + Bk+ig'l+i.k' + Bk'+iglk'+i + Pk-ig'Li.k- + Pk'-ig'L'-i) 1^ + ^Kk'. (V.iOa) 
5fe,fe' = {QkM'g'k.k' - Bk+ig'k+i,k' - Bk'+ig'k.k'+i - Pk-ig'k-i.k' - Pk'-ig'k.k'-i) (V.iOb) 

(for fc, fc' > 2) 

In the high frequency limit, (jV.lOp reduces to 

^ oo) = i^fc^fc', (V.lla) 

gl,,{iu^^)^0. (V.llb) 

The shear stress ak^k' ascribed to the (fc, fc')-chain is derived from the following formula for the Gaussian chain: 

CTfc.fc' -^TT f xyFk^k'ir,t) = e [G'k^i,,{uj) sin ujt + G'lj,,{uj) cos cut] (V.12) 



where the dynamic shear moduli with respect to the (fc, fc')-chain are defined by 

Gl,,{iu)^kBTgik'{^), (V.13a) 
Gl,,{iu) ^ kBTgl^,{uj). (V.13b) 

Since we are assuming that any chains whose both ends are associated with other functional groups are elastically 
effective, the total moduli are given by 



G'(^)=fcBT^^4,,(c.), (V.14a) 

fe>2 k'>2 

G"{u:)=kBTY,Y.9U{^)- (V.14b) 



k>2 k'>2 



^ Ea.l lV.9l l assures that the number of (fc, fc')-chains of the first order, say u^^y{t), is zero; '^l^l,/{t) = f '^^^k^ii^^'t) = 0- 
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2P 2P 

FIG. 10: The breakage rate of a pairwise junction is 2/3 because it decays if one of two functional groups in the junction is 
disconnected from the other. Since active chains are connected with two pairwise junctions through both ends, the breakage 
rate of the active chain is 4/3. 



The high-frequency plateau modulus is then found from (jV.lla[) to be 

Goo =G'{lu^oo)^ lyfksT {V.15) 

where 

k>2 k'>2 

is the total number of active chains. 

As an example, let us study the simplest case that the association is allowed only in pairs. In other words, only 
junctions whose multiplicity is two are allowed to exist in the system. There is no difference between two models for 
junctions (i.e., the saturating junction model and fixed multiplicity model) in this case. Since we are treating telechelic 
polymers, only linearly extended chains can be formed rather than a three dimensional network. Nevertheless, it is 
worth studying such simple situation because this is the only case that the moduli can be analytically expressed in 
simple forms, thereby giving insight into the system with multiple junctions. Note that 53 2 = ff2 3 = 2 = 52 3 — 
since there are no 3-j unctions. Then (jV.lOp reduces to 

52,2 = -Q^.ig'i.il'^ + ^2,2, (V.17a) 
52,2 = Qia9i,2l<^- (V.17b) 
As a result, the moduli becomes the Maxwellian with a single relaxation time 1/Q2,2, i.e., 



i2 



G'/ksT ^ .9^.2 = ' (V-18a) 

\Q2a^) + -1- 

C'lksT = g'i, = ^i^^z^f , (V.18b) 
{Q2.2^Y + 1 

where the number of active chains is v^^ = 1^2,2 = n{zhiqi)^ , and the reciprocal of the relaxation time, or the breakage 
rate of the active chain, is given by Q2,2 = 4/3-I-2P2 = 4/3.' The coefficient 4 of the chain breakage rate stems from the 
annihilation rate 2(3 of the junction to which one end of the chain is belonging and 2(3 for the other end (see Fig fTU)) . 



A. Saturating Junction Model 

Now we study dynamic- mechanical and viscoelastic properties of the system where multiple junctions (i.e., k > 2) are 
allowed to exist. We first consider the saturating junction model, i.e., the junction takes a limited range k — 1,2, ...Sm 
of the multiplicity. This condition can be attained by employing (|IV.10p for P^., i.e., 

_ ( (3zk (1 < k < Sm " I) (.r -, Q^, 

^'^-\o (fc>.„) ■ ^^■^^> 

By putting (jV.lQP into (jV.lOp and by solving a set of linear algebraic equations for g'f. j., and g'l j^, , we can obtain the 
dynamic shear moduli with the help of (|V.14p where the summation is taken over 2 < fc, fc' < s™. Note that the number 
of unknowns in (jV.10[l is Sm{sm ~ 1); when — 4, for example, there are 12 unknowns: g'2 2j 53 2' 53 31 54 2j 54 31 54 4 

and 32,2: 53,2; 53,3! 54,2: 54,3' 54,4- 



The probability p2 that an unroacted group is connected with a 2-junction is zero since 3-junction is not allowed to exist by definition. 
Therefore, P2 oc p2 = 0. 
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c c 

FIG. 11: (i) The dynamic shear moduh (reduced by nksT) for the saturating junction model as a function of the frequency. The 
reduced polymer concentration c is varying from curve to curve, while the maximum multiplicity is fixed at Sm = 10. (ii) The 
reduced plateau modulus, (iii) reduced zero-shear viscosity, and (iv) relaxation time plotted against the reduced concentration. 
The maximum multiplicity is varying from curve to curve. The insets of (ii) ~ (iv) show the linear-log plot of each quantity as 
a function of c. 




^ increasing c 

FIG. 12: The breakage rate of the chain linking two k = 2 junctions (upper figure) is 4/3 (see Fig llOp . On the other hand, the 
breakage rate of the chain connecting two junctions with the multiplicity fc > 3 (lower figure) is 2/3 because the chain debridges 
if its one of two ends are dissociated from the junction. 



Fig[TT] (i) shows the dynamic shear moduli reduced by nksT as a function of the unitless frequency. The reduced 
polymer concentration is changed from curve to curve for a maximum multiplicity fixed at Sm = 10. It appears that 
they are Maxwellian with a single relaxation time. To see the detail, the plateau modulus Goo — linia;^oo G" (w) 
reduced by nksT is drawn in Fig llll (ii) as a function of c. The reduced plateau modulus increases as irrespective 
of the value of Sm when c is small™ and approaches to 1 with increasing c. Fig llll (iv) shows the relaxation time 
determined from the peak position of G" as a function of c. The relaxation time becomes close to 1/ (4/3) in the limit of 
low reduced concentration because most active chains belong to junctions with A: = 2 in this limit. On the other hand, 
T approaches to 1/(2/3) with increasing c because more active chains tend to connect with junctions with fc > 3 when 
c is large (see Fig |12l) . Thus t increases with increasing c. The zero-shear viscosity rjo = limj^^o G"{uj)/uj reduced by 
nksT/ (3 is shown in Fig llll (iii) as a function of c. The zero-shear viscosity is roughly estimated to be 770 ~ GooT (or 
rjo/ (nksT/ (3) ~ Goo/(?^fcBT)/3r), and hence the reduced viscosity starts to increase at c = and approaches to 0.5 
with increasing c. 

Fig |13l (i) shows the reduced shear moduli for several maximum multiplicity s,„. The reduced plateau modulus. 



When c is small, the extent of reaction is proportional to c, i.e., « ~ c/(l + c) ~ c. Therefore, we find Gaa /^kgT = ~ c^. 
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FIG. 13: (i) The dynamic shear moduh (reduced by nfcsT) for the saturating junction model as a function of the frequency. 
The maximum multiphcity Sm is varying from curve to curve for a reduced concentration fixed at c = 5. (ii) The reduced 
plateau modulus, (iii) reduced zero-shear viscosity, and (iv) relaxation time plotted against the mgiximum multiplicity for 
several reduced concentrations. 



reduced zero-shear viscosity and relaxation time are also plotted in Fig |13l as a function of s„j. Let us comment 
on the relaxation time here. The relaxation time is determined from the ratio between the number of active chains 
incorporated into the junction with multiplicity fc > 3 and that with k = 2. When Sm — 2, for example, all active 
chains belong to 2-junctions, and therefore we have t = 1/(4/3) irrespective of the value of c. When Sm > 3, more 
chains can connect to junctions with A: > 3, so that the fraction q2 of chains belonging to 2-junctions decreases (this 
tendency can also be confirmed from FiglS]) . Thus r increases with increasing Sm and approaches to a fixed value for 
each c. As for the dependence of Goo/nksT = o? on Sm, see IIV Al for reference. 



B. Fixed Multiplicity Model 

Next, we consider the fixed multiplicity model, i.e., the multiplicity can take only = 1 (unreacted) and k = s 
(reacted) for all junctions. This condition can be approximately attained by employing (|IV.20p for P^, i.e., 

( S (1 < A: < s - 1) 

Pfc = pzk \ (5-(^-2) (fc = .s - 1) , (V.20) 
[ (fc > s - 1) 

where we are setting S = 0.01. By putting (jV.20[) into (jV.10[) and by solving a simultaneous equation, we can obtain 
the dynamic shear moduli for the (approximate) fixed multiplicity model. 

Fig 1 141 shows the dependence of the reduced shear moduli on c. The reduced plateau modulus (figure (ii)) increases 
with increasing s for c > 2, but it decreases for c < 2 according to the s dependence of a (see FiglS]). The relaxation 
time (figure (iv)) is almost 1/(2/3) irrespective of the values of s(> 3) and c because almost all junctions take the 
multiplicity larger than or equal to three. However, there is a small number of 2-junctions in the system, especially 
when c is close to 0, due to the finite 6. Thus, the relaxation time becomes smaller than 1/(2/3) (but larger than 
1/(4/3)) for c close to 0. The reduce zero-shear viscosity (figure (iii)) is approximately half the reduced plateau modulus 
because /3t is about 0.5. 
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FIG. 14: (i) The dynamic shear moduli (reduced by nksT) for the fixed multiplicity model as a function of the frequency. The 
reduced concentration is varying from curve to curve while the multiplicity is fixed at s = 5. (ii) the reduced plateau modulus, 
(iii) reduced zero-shear viscosity, and (iv) relaxation time plotted against the reduced concentration. 



VI. SUMMARY AND DISCUSSIONS 



In this paper, we developed a theory for transient networks with muhiple junctions. One of our objectives is to 
investigate relations between rheological (and equilibrium) properties of the network and the junction multiplicity 
in the absence of looped chains. We assumed that the connection rate of a functional group is proportional to the 
volume fraction of junctions to which it is going to connect, and showed that the law of mass action holds. As an first 
attempt, we defined active chains locally, i.e., chains whose both ends are connected with any other functional groups 
are elastically effective. The dynamic shear moduli are well described in terms of the Maxwell model characterized by 
a single relaxation time and the high-frequency plateau modulus (and the zero-shear viscosity). The reduced moduli 
depend on thermodynamic quantities such as polymer concentration and temperature through the reduced polymer 
concentration c. The plateau modulus and the zero-shear viscosity increase nonlinearly at small c and proportional 
to c when c is large. The relaxation time also increases with increasing c due to the presence of pairwise junctions at 
small c. 

The modulus, viscosity and relaxation time start to rise at c = 0. This indicates that the critical concentration for 
the sol/gel transition is c* = in this model. In other words, the system is always in the postgel regime irrespective 
of the amount of polymers. This unfavorable result can be ascribed to the lack of global information in the current 
definition of active chains. 

The present theory reduces to the TE theory in the high reduced concentration limit c ^ oo. Let us take a 
summation over 2 < k,k' < s,n in pi.7|) . Then we have 



+ V • (rF(r, t)) = -2/3(r)F(r, t) - B^it) ^ (F,,2(r, t) + F,,kir, t)) 

fe=2 

Sm-1 

+(E (vi.i) 

k=l 

where 

F(r,i)^^^Ffc.fc,(r,i) (VL2) 



fc=2 fc'=2 
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logc logc 

FIG. 15: (i) Effects of the terms in (|VI.1|) and (|VI.4|) associated with the creation/annihilation process of 2-junctions, and (ii) 
the fraction of chains belonging to 2-junctions plotted against the reduced polymer concentration. The maximum multiplicity 
is varying from curve to curve in both figures. 



is the total number of active chains with the head-to-vector r, and 

\t) = Y,{vkAt) + vi.k{'t)) (VI.3) 



k=2 



is the total number of dangling chains. (We are assuming that the dissociation rate does not depend on the junction 
multiplicity as in the text.) The second term of the right-hand side in (jVI.l[) and pi inside the second parentheses 
are related with the annihilation/creation process of 2-junctions. As shown below, the relative amount of 2-junctions 
becomes negligible in the high c limit, and hence effects of these terms disappear. Let us consider the case that a small 
shear deformation is applied to the system as discussed in the text. Upon integration with respect to r, eg. ljVI.ip 
becomes 

= -2p{iy'ff + ^ka) + iY.P^+ Pi) ■ 



fc>2 k>l 

The terms X]fc>2 ^k,2 and pi associated with 2-junctions satisfy the following relation: 

Yjk>2 ^k,2 pi q2 



J2k>lPk Oi 



(VI.5) 



With increasing the reduced concentration, (?2/a approaches to zero as shown in Fig[T5](i)," and hence the right-hand 
side of (|VI.4p becomes close to —2Pu'^f^ + J2k>iPk- K indicates that (|VI.ip reduces to 

+ V • (rF(r, t)) = -2pF{r, t) + pv'' h(v) (VI.6) 

in the high c limit where 

P = ^Pfe (VI.7) 

k>l 

is the probability that an unreacted group connects to any junctions per unit time. Eq. (jVf .6^ is equivalent to the 
fundamental equation of the TE theory when we regard p as a constant and assume that isolated chains are absent." 

In Fig ll6l theoretically obtained plateau modulus, zero-shear viscosity and relaxation time for the saturating junc- 
tion model are compared with with experimental data for aqueous solutions of telechclic PEO end-capped with Cig 
alkanes (sl, The reduced concentration c used in the theory was converted into the polymer concentration in 
weight percentage Cw through a relation c = ^c^,, where ^ = (2000iVyi/M)Awo {Na is Avogadro's number). We see 



" Eq. ifVLSt can be approximately expressed as 1/(1 c) when Sm 

is large. 

° 2/3 is the transition rate from active chains to dangling chains. This quantity is denoted as /3 in the TE theory l2Cl l2ll . 



22 



(i) 



(ii) 




5000 




polymer concentration [wt%] 

FIG. 16: Comparison between theoretical results for the saturating junction model and experimental data obtained for (i) 
telechelic PEO with narrow molecular weight distributions {Mw = 35kg/mol) fully end-capped with Cie alkanes reported by 
Pham et al. 9] (called HDU-35) and (ii) hydrophobically modified ethylene oxide-urethane copolymers (HEUR) with similar 
molecular weight end-capped with the same hydrophobes reported by Annable et al. 0) (called hd-35 after ref.[25l|). Theoretical 
curves for the zero-shear viscosity and the relaxation time are drawn with /3 = 90 (1/sec) for IfDU-35 and (5 = 3.5 (1/sec) for 
hd-35Q A factor ^ (see text) is assumed to be 1 for HDU-35 and 0.35 for hd-35, respectively. 

"A discrepancy in the value of /3 for the same alkanes might stem from the difference in the coupling agents between the alkanes and the 
PEO backbone [llil. 



that both agree fairy weU with each other for larger Sm except for the relaxation time of HEUR aqueous solutions. 
According to Annable et at. , a concentration dependence of the relaxation time is strongly affected by the amount 
of superbridges formed by connecting several primary bridges in series because they have shorter lifetime due to a 
number of possible disengagement points inside. Actually, we can infer from Fig |15l (i) that a considerable fraction 
of active chains are forming superbridges at low concentration; a fraction of active chains ((72/0:)^ whose both ends 
are connected with 2-junctions becomes close 1 when c — » although both a and (72 approach to (see FiglS] (ii) 
and Fig llSI (ii) respectively). In order to treat such superbridges in more appropriate manner, we need to consider a 
global structure of the network by introducing a concept of the path connectivity to the infinite network according to 
refs. [S^, [35| . The effect of superbridges on the viscoelasticities will be discussed in detail in the following paper. 



APPENDIX A: DERIVATION OF (|II.30|) 



The number of the state X depicted in FigllT] (corresponding to the state A in FigOwith kh — 1 and kt — 0) is 

^2,fc'(r,rW;i) = ^F,.k,{r,t)f^''\r^'^\t). (A.l) 
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FIG. 17: Transitions between the k = 1 chain and the k — 2 chain. 



The (1, A;')-chain is generated from the state X if the head of the (2, fc')-chain or the head of another chain is detached 
from the 2-junction (see the transition from the state X to state Y in Fig fTT]) . Then the increment of (1, fc')-chains 
due to this process is estimated to be 



dr (/32 (r) + /32 (r J-2,fe' (r, r^'^) ; t) 



(A.2) 



Similarly, the increment of (1, A;')-chains as a result of the transition from the state X' (in FigfTT]) to Y' is written as 



(A.3) 



2 dr f32{r)F2,k'ir,t) + - ( j dr I32{r)f^"{r,t) ] v^^it)- 



By adding (|X2|) and ([XS]) . we obtain (|lL30l) . 
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